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ABSTRACT 



Aims. A transient astrophysical event observed in both gravitational wave (GW) and electromagnetic (EM) channels would yield 
rich scientific rewards. A first program initiating EM follow-ups to possible transient GW events has been developed and exercised 
by the LIGO and Virgo community in association with several partners. In this paper, we describe and evaluate the methods used to 
promptly identify and localize GW event candidates and to request images of targeted sky locations. 

Methods. During two observing periods (Dec 17 2009 to Jan 8 2010 and Sep 2 to Oct 20 2010), a low-latency analysis pipeline was 
used to identify GW event candidates and to reconstruct maps of possible sky locations. A catalog of nearby galaxies and Milky Way 
globular clusters was used to select the most promising sky positions to be imaged, and this directional information was delivered to 
EM observatories with time lags of about thirty minutes. A Monte Carlo simulation has been used to evaluate the low-latency GW 
pipeline's ability to reconstruct source positions correctly. 

Results. For signals near the detection threshold, our low-latency algorithms often localized simulated GW burst signals to tens 
of square degrees, while neutron star/neutron star inspirals and neutron star/black hole inspirals were localized to a few hundred 
square degrees. Localization precision improves for moderately stronger signals. The correct sky location of signals well above 
threshold and originating from nearby galaxies may be observed with ~5Q% or better probability with a few pointings of wide-field 
telescopes. 



Key words, gravitational waves - methods: observational 
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1. Introduction 

The Laser Interferometer Gravitational- Wave Obser vatory 
(LIGO) (Abbott etal. 2009a) and Virgo dAccadia et all I20T1I) 
have taken significant steps toward gravitational wave (GW) 
astronomy over the past decade. The LIGO Scientific 
Collaboration operates two LI GO observatories in the U.S. 
along with the GEO 600 detector dLiick et al.ll2010l) in Germany. 
Together with Virgo, located in Italy, they form a detector net- 
work capable of detecting GW signals arriving from all direc- 
tions. Their most recent joint data taking run was between July 
2009 and October 2010. GEO 600 and Virgo are currendy oper- 
ating during summer 201 1, while the LIGO interferometers have 
been decommissioned for t he upgrade to the next-generation 
Advanced LIGO detectors dHarrv et al.l 1201 Oh . expected to be 
operational around 2015. Virg o will also b e upgraded to be- 
come Advanced Virgo ( Acerne se et al. Il2008h. Additionally, the 
new LCGT detector dKuroda & The LCGT Collaborationll2010l) 
is planned in Japan. These "advanced era" detectors are expected 
to detect compact binary coalescences, possibly at a rate of 
dozens per year, after reaching design sensitivity dAbadie et alJ 
l20To3) . 

Detectable, transient GW signals in the LIGO/Virgo fre- 
quency band require bulk motion of mass on short time scales. 
Emission in other channels is also possible in many such rapidly 
changing massive systems. This leads to the prospect that some 
transient GW sources may have corresponding electromagnetic 
(EM) counterparts which c ould be discover e d with a low latenc 
response to GW triggers ( Svlvestre 2003 ; Kanner et al. 200 
IStubbsll2008l:lKulkarni & Kasliwalll2009l: iBloom et al .I l2009h . 

Finding these EM counterparts would yield rich scientific 
rewards (see Sect. O, but is technically challenging due to im- 
perfect localization of the gravitational wave signal and uncer- 
tainty regarding the relative timing of the GW and EM emis- 
sions. This paper details our recent effort to construct a prompt 
search for joint GW/EM sources between the LIGO/Virgo de- 
tector network and partner EM observatories (see Sect. [3]). The 
search was demonstrated during two periods of live LIGO/Virgo 
running: the "winter" observing period in December 2009 and 
January 2010 and the "autumn" observing period in September 
and October 2010. We focus here on the design and performance 
of software developed for rapid EM follow-ups of GW candi- 
date events, as well as the procedures used to identify significant 
GW triggers and to communicate the most likely sky locations 
to partner EM observatories. The analysis of the observational 
data is in progress, and will be the subject of future publications. 



2. Motivation 

2.1. Sources 

A variety of EM emission mechanisms, both observed and the- 
oretical, may occur in association with observable GW sources. 
Characteristics of a few scenarios helped inform the design and 
execution of this search. Here, some likely models are presented, 
along with characteristics of the associated EM emission. 

2.LL Compact Binary Coalescence 

Compact binary systems consisting of neutron stars and/or black 
holes are thought to be the most common sources of GW emis- 
sion detectable with ground-based interferometers. Radiation of 
energy and angular momentu m causes the orbit to decay (in- 
spiral) until the objects merge dCutler et al.lll993l ). For a system 



consisting of two neutron stars (NS-NS) or a neutron star and 
a stellar-mass black hole (NS-BH), the inspiral stage produces 
the most readily detectable GW signal. The energy flux reach- 
ing Earth depends on the inclination angle of the binary orbit 
relative to the Une of sight. The initial LIGO- Virgo network is 
sensitive to optimally oriented NS-NS mergers from as far away 
as 30 Mpc, an d mergers between a NS and a 10 Mq black hole 
out to 70 Mpc dAbadie et alj|2010d) . Models of the steflar com- 
pact object population in the local universe estimate the rate of 
NS-NS mergers detectable with initial detectors to be between 
2 X 10""* and 0.2 per year With advanced detectors, these range 
limits are expected to increase to 440 and 930 Mpc, respectively. 

The energetics of these systems suggest that an EM counter- 
part is likely. The final plunge radiates of order lO^^* ergs of grav- 
itational binding energy as gravitational waves. If even a small 
fraction of this energy escapes as photons in the observing band, 
the resulting counterpart could be observable to large distances. 
The EM transients that are likely to foUow a NS-BH or NS-NS 
merger are described below. 

Short-hard gamma-ray bursts (SGRBs), which typically have 
durations of 2 seco nds or less, may be powered by NS-NS 
or NS-BH mergers dNakaijl2007l: iMeszaros 2006: Piran2004 . 
Afterglows of SGRBs have been observed in wavelengths 
from radio to X-ra y, and out to Gpc distances dNakad 120071 : 
iGehrels et alJ l2009l) . Optical afterglows have been observed 
from a few tens o f seconds to a few d ays after the GRB trigger 
(see, for example, iKlotz et al] d2009al) ). and fade with a power 
law r", where a is between 1 and 1.5. At 1 day after the trigger 
time, the apparent optical magnitude would be between 12 and 
20 for a source at 50 Mpc (Kann et al. 201 1), comparable to the 
distance to which LIGO and Virgo could detect the merger 

Even if a compact binary coalescence is not observable in 
gamma-rays, there is rea son to expect it will pro duce an observ- 
able optical counterpart. iLi & Paczvhskil d 19981) suggested that, 
during a NS-NS or NS-BH merger, some of the neutron star's 
mass is tidally ejected. In their model, the ejected neutron-rich 
matter produces heavy elements through r-process nucleosyn- 
thesis, which subsequently decay and heat the ejecta, powering 
an optical afterglow known as a kilonova. The predicted optical 
emission is roughly isotropic, and so is observable regardless of 
the orientation of the original binary system. This emission is 
expected to peak af ter about one day, ar ound magnitude 18 for 
a source at 50 Mpc dMetzger et alj|20l"ol) . and then fade over the 
course of a few days following the merger. 

2.L2. Stellar Core Collapse 

Beyond the compact object mergers described above, some other 
astrophysical processes are plausible sources of observable GW 
emission. GW transients with unknown waveforms may be dis- 
covered by searching the LIGO and Virgo data for short periods 
of excess power (bursts). The EM counterparts to some likely 
sources of GW burst signals are described below. 

Core-collapse supernovae are likely to produce some amount 
of gravitational radiation, though large uncertainties still exist 
in the expected waveforms and energetics. Most models predict 
GW spectra that would be observable by initial LIGO and Virgo 
from distances within some fraction of the Milky Way, but not 
from th e Mpc di stances needed to observe GWs from another 
galaxy dOtll2009l) . Neutrino detectors such as SuperKamiokande 
and IceCube should al so detect a large number of neutrinos from 
a Galactic supernov a ( Beacom & Vogeliri999l : iHalzen & Raffeltl 



120091: ILeonor et al.ll2010h . Galactic supernovae normally would 
be very bright in the optical band, but could be less than obvious 
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if obscured by dust or behind the Galactic center. Optical emis- 
sion would first appear hours after the GW and neutrino signal 
and would peak days to weeks later, fading over the course of 
weeks or months. 

Long-soft gamma-ray bursts (LGRBs) are b elieved to be as- 
sociated with the core collapse of massive stars ('Wooslev''l993t 
MacFadyej, &_ Wgogley 1 999; Piran 2004; Woosley & Bloonj 



dCutler et al.lll993l;IVanisnerill2000l;lFlanagan & Hindereiil2008l; 



Andersson et al. ~ |201 U iPannarale et al.l 1201 U iHinderer et al.l 
2010h . 



2006t iMetzger et ai.l 120111) . A large variety of possible GW 



emitting mechanisms within these systems have been pro- 
posed, with some models predicting GW spectra that would 



and Virso dFrver et al. 
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Corsi & MeszarosI 


2009 
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). The afterglows of LGRBs, like the 



afterglows of SGRBs, typically show power law fading with 
a = 1 - 1.5. However, the peak isotropic equivalent luminos- 
ity of LGRB after glows is typically a factor of 1 brighter than 
SGRB afterglows dN^l20()7HKann et alllloTol) . 

An off-axis or sub-energetic LG RB ma y also be observed as 
an orp han afterglow or dirty fireball tGranot et a l. 2002; Rho adsl 
12001 . These transients brighten over the course of several days 
or even weeks, depending on the observing band and viewing an- 
gle. Identifying orp han afterglows in large area surveys, such as 
iRvkoffet al.l(l2005h . has proven difficult, but a GW trigger may 
help distinguish orphan afterglows from other EM variability. 



2.1.3. Other Possible Sources 

Cos mic string cusps are another possible joint s ource of 
GW ("Siemens et alJ 120061: iDamour & VilenkinI l2000l) and EM 
jVac haspati 2008) radiation. If present, their distinct GW sig- 
nature will distinguish them from other sources. On the other 
hand, even unmodeled GW emissions can be detected using GW 
burst search algorithms, and such events may in some cases pro- 
duce EM radiation either through internal dynamics or through 
interaction with the surrounding medium. Thus, our joint search 
methods should allow for a wide range of possible sources. 



2.2. Investigations enabled by joint GW/EM observations 

A variety of astrophysical information could potentially be ex- 
tracted from a joint GW/EM signal. In understanding the pro- 
genitor physics, the EM and GW signals are essentially com- 
plementary. The GW time series directly traces the bulk motion 
of mass in the source, whereas EM emissions arising from out- 
flows or their interaction with the interstellar medium give only 
indirect information requiring inference and modeling. On the 
other hand, observing an EM counterpart to a GW signal reduces 
the uncertainty in the source position from degrees to arcsec- 
onds. This precise directional information can lead to identifica- 
tion of a host galaxy, and a measurement of redshift. Some spe- 
cific questions that may be addressed with a collection of joint 
GW/EM signals are discussed below. 

If the GW source is identified as a NS-NS or NS-BH merger, 
additional investigations are enabled with an EM counterpart. 
The observation of the EM signal will improve the estimation 
of astrophysical source parameters. For example, when attempt- 
ing parameter estimation with a bank of templates and a sin- 
gle data stream, the source's distance, inclination angle, and 
angular position are largely degenerate. A precise source po- 
sitio n from an EM counterpart would help break this degener- 
acy (iDalal et al.ll2006HNissanke et al.ll2010l) . High precision pa- 
rameter estimation may even constrain the NS equation of state 



Observing EM counterparts of NS-NS and NS-BH merger 
events will give strong evidenc e as to which class of source, if 
either, is the source of SGRBs dBloom et al ] l20()9h . In addition, 
if some neutron star mergers are the sources of SGRBs, a collec- 
tion of joint EM/GW observations would allow an estimate of 
the SGRB jet opening angle by comparing the number of merger 
events with and without observable prompt EM emission, and 
some information would be obtainable even from a single loud 
event (Kobayashi & Meszaros 2003b; Seto 200"^). 

An ensemble of these observations could provide a novel 
measurement of cosmological parameters. Analysis of the well- 
modeled GW signal will provide a measurement of the lumi- 
nosity distance to the source, while the redshift distance is mea- 
surable from the EM data. Taken together, they provide a di- 
rect measurement of the local Hubble constant (ISchutzlll986 l : 
iMarkoviJ 1991 IDalal et al]l2006l:lNissanke et al.ll201flh 

Finally, all of the above assume that general relativity is the 
correct theory of gravity on macroscopic scales. Joint EM/GW 
observations can also be used to test certain predictions of gen- 
eral relativ ity, such as the propagation speed and p olarizations of 
GWs (IWillll2005l:lYunes et alJl20Tol:lKahvall201 Ih . 

In the case that the transient GW source is not a binary 
merger event, the combination of GW and EM information will 
again prove very valuable. In this scenario, the gravitational 
waveform will not be known a priori. Any distance estimate 
would be derived from the EM data, which would then set the 
overall scale for the energy released as GWs. 

As in the merger case, the linking of a GW signal with a 
known EM phenomenon will provide insight into the underly- 
ing physical process. For example, the details of the central en- 
gine that drives LGRBs are unknown. The GW signal could give 
crucial clues to the motion of matter in the source, and poten- 
tially distinguish between competing models. A similar insight 
into the source mechanism could be achieved for an observa- 
tion of GW emission associated with a supernova. Rapid identi- 
fication may also allow observation of a supernova in its ear- 
liest moments, an opp ortunity that currently depends on luck 
jSoderberg et al.ll2008l) . 

2.3. Extend GW Detector Reach 

Finding an EM counterpart associated with a LIGO/Virgo trigger 
would increase confidence that a truly astrophysical event had 
been observed in the GW data. Using EM transients to help dis- 
tinguish low amplitude GW signals from noise events allows a 
lo wering of the detection threshold, as was done in searches such 
as I Abbott et"aD (l2010l) . iKochanek & PiranI (1 19931) estimated that 
the detectable amplitude could be reduced by as much as a fac- 
tor of 1.5, increasing the efifective detector horizon distance (the 
maximum distance at which an optimally oriented and located 
system could be detected) by the same factor and thus increas- 
ing the detection rate by a factor of 3. In practice, the actual im- 
provement in GW sensitivity achieved by pairing EM and GW 
observations depends on many factors unique to each search, in- 
cluding details of the source model and data set, and so is diffi- 
cult to predict in advance. 

In the case of a coincidence between a GW signal and a dis- 
covered EM transient, the joint significance may be calculated 
by assuming that the backgrounds of the EM and GW search are 
independent. The False Alarm Rate (FAR) of a GW/EM coin- 
cidence is the FAR of the GW signal, as described in Sect. 14.21 
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times a, the expected fraction of observations associated with a 
false or unrelated EM transient. The false alarm fraction a may 
be estimated using fields from surveys not associated with GW 
triggers. The measured value of a will depend heavily on the 
telescope being used, the cuts selected in image analysis, galac- 
tic latitude of the source and other factors. For example, surveys 
with the Palomar Transient Factory require a sophisticated clas- 
sification mechanism for rejecting contaminants. Each set of im- 
age subtractions covering 100 - 200 square degrees yields ~ 10^ 
candidates. Of these, 30 - 150 sources are selected after impos- 
ing cuts optimized for the detection of fast evolving transients 
(Bloom et al. 201 1]) . Using classificati on software designed for 
PTE data (Oarical) (iBloom et al.ll20Tlh . the selected sources un- 
dergo an automatic classification as type "transient" or "variable 
star" based on time-domain and context properties. Promising 
candidates are selected for additional, spectroscopic observa- 
tion. Of the sources that are classified as transients, and then fol- 
lowed up spectroscopically, ~ 82% are supernovae dBloom et alJ 
l20Tlh . To use EM transience to improve confidence in a GW 
signal, the time-domain sky in the wavelength of interest must 
be well understood. Transients that are found in directional and 
time coincidence with GW triggers would increase confidence 
in the GW signal only if the cha nce of a similar, incidental c oin- 
cidence is understood to be low dKulkarni & Kasliwalll2009h . 

2.4. Implications for Search Design 

Characteristics of the target sources helped determine when and 
where to seek the EM counterparts to GW event candidates. Eor 
reasons discussed in this section, the search strategy presented 
in this paper emphasizes capturing images as soon as possible 
after the GW trigger, along with follow-up images over subse- 
quent nights. The rates of stellar core-collapse and compact ob- 
ject mergers within our own galaxy are much less than one per 
year, and so field selection was strongly weighted towards re- 
gions containing nearby galaxies. 

The observations and theoretical models of EM transients 
discussed above provided guidance when choosing the observ- 
ing cadence. GRB optical afterg lows have been ob served dur- 
ing the prompt emission phase ( lKlotzetalJl2009bl) and up to 
many hours after the trigger. Eor this search, the first attempt to 
image the source position was made as soon as possible after 
validating a GW trigger. In both the kilonova dLi & Paczvriskil 
Il998h and supernova ( Ott 2009) models, some time lag exists 
between the release of GW and EM emission, primarily due to 
the time it takes the outflowing material to become optically thin. 
This time lag may be from several hours for a k ilonova, up to 
days f or a core-collapse supernova. Furthermore. 'Cow ard et al.l 
(1201 lb show that for GRBs that are off-axis, the optical afterglow 
may not be visible until days after the burst. For these reasons, 
repeated observations over several nights are desirable. Also, the 
light curves obtained by observing the same fields over multiple 
nights are critical clues for discovering and classifying transient 
sources. 

Knowing where to look for the counterpart to a GW trig- 
ger is challenging. Directional estimates of low signal to noise 
ratio (SNR) binary inspiral sources with the 2009-10 GW detec- 
tor network hav e uncertainties of several tens of square degrees 
(lFairhursdl2009l) . This suggests using telescopes with a field of 
view (FOV) of at least a few square degrees if possible. Even 
with such a "wide field" instrument, there is a striking mismatch 
between the large area needing to be searched, and the size of a 
single FOV. 



The problem may be partially mitigated by making use of 
the known mass distribution in the nearby universe. A search for 
GW counterparts can dramatically reduce the needed sky cov- 
erage by focusing observa tions on galaxies within the distance 
limits of the GW detectors (iKanner et al.ll2008l:lNuttall & SuttonI 
120101) . Limiting the search area to known galaxies may also 
improve the feasibility of identifying the true counterpart 
from among other objects with time-varying EM emissions 
(iKulkarni & KasUwall [20091) . Even within the Milky Way, a 
search may emphasize known targets by seeking counterparts 
within globular clusters, where binary systems of compact ob- 
jects may form efliciently (O'Learv et al. 2007). 

An emphasis on extragalactic and globular-cluster sources 
has the potential drawback that any counterparts in the plane of 
the Milky Way may be missed. Also, neutron star mergers that 
occur at large distances from their host galaxies may not be ob- 
served, thou^hJh£_go£u^ large kicks should be small 
dBergedl20ldlKeUev et al.ll2010l) . 

Our selection of fields to observe was weighted towards ar- 
eas containing known galaxies within 50 Mpc. The utilized cata- 
log of nearby galaxies and globular clusters, and the process for 
selecting fields to observe, is described in Sect. |5] 

3. GW and EM Instruments 

3.1. Gravitational Wave Detector Network 

The LIGO and Virgo detectors are based on Michelson-type 
interferometers, with Fabry-Perot cavities in each arm and a 
power recycling mirror between the laser and beamsplitter to 
dramatically increase the power in the arms relative to a simple 
Michelson design. The GEO 600 detector uses a folded interfer- 
ometer without Fabry-Perot arm cavities but with an additional 
recycling mirror at the output to resonantly enhance the GW sig- 
nal. As a gravitational wave passes through each interferometer, 
it induces a "strain" (a minuscule change in length on the order 
of 1 part in 10^' or less) on each arm of the interferometer due 
to the quadrupolar perturbation of the spacetime metric. The in- 
terferometers are designed to measure the differential strain on 
the two arms through interference of the laser light when the 
two beams are recombined at the beam splitter, with the rela- 
tive optical phase m odulated by the passing gravitational wave 
(lAbbott et al.ll2009ah . 

In 2009-2010 there were two operating LIGO interfer- 
ometers, each with 4-km arms: HI, located near Hanford, 
Washington, and LI, located in Livingston Parish, Louisiana^ 
Virgo (VI) has arms of length 3 km and is located near Cascina, 
Italy. GEO 600 data was not used in the online search described 
in this paper, but was available for offline reanalysis of promis- 
ing event candidates. The large physical separation between the 
instruments means that the effects of local environmental back- 
ground can be mitigated by requiring a coincident signal in mul- 
tiple interferometers. Each interferometer is most sensitive to 
GW signals traveling parallel or anti-parallel to zenith, but the 
antenna pattern varies gradually over the sky, so that the detec- 
tors are essentially all-sky monitors. 

The EM follow-up program described in this paper was exer- 
cised during the 2009-2010 science runs. While single-detector 
triggers had been generated with low latency in earlier science 

' Earlier science runs included a second interferometer at Hanford, 
called H2, with 2-km arms. H2 will reappear as part of Advanced LIGO, 
either as a second 4-km interferometer at Hanford or else at a site in 
Western Australia. The latter option w ould greatly improve the sour ce 
localization capabilities of the network dFairhursteOl lh[Schutzll201 Ih . 
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runs for diagnostic and prototyping purposes, 2009-2010 was 
the first time that a systematic search for GW transients using 
the full LIGO- Virgo network was performed with low latency, 
and the first time that alerts were sent to external observatories. 



3.2. Optical and Other Electromagnetic Observatories 

In an effort to explore various approaches, the telescope network 
used in 2009-10 was intentionally heterogeneous. However, 
most of instruments had large fields of view to accommodate 
the imprecise GW position estimate. The approximate location 
of each EM observatory is shown in Fig.[T] and Tables [1] and |2] 
show some of the properties of each observatory. 



3.2.1. Optical Instruments 



The Palomar Transient Factory (PTF) jLaw et al. ] l2009t 
iRau et all l2009h operates a 7.3 square degree FOV camera 
mounted on the 1.2 m Oschin Telescope at Palomar Observatory. 
A typical 60 s exposure detects objects with a limiting magni- 
tude R = 20.5. For the autumn observing period, the PTF team 
devoted ten fields over several nights at a target rate of 1 trigger 
for every three weeks. 

Pi of the Sky (iMalek et al.ll2009l) observed using a camera 
with a 400 square degree FOV and exposures to limiting mag- 
nitude 11-12. It was located in Koczargi Stare, near Warsaw. 
The camera was a prototype for a planned system that will si- 
multaneously image two steradians of sky. The target rate was 
approximately 1 per week in the autumn run, followed up with 
hundreds of 10 s exposures over several nights. 

The QUEST camera (iBaltav et al.ll2007h . currently mounted 
on the 1 m ESO Schmidt Telescope at La Silla Observatory, 
views 9.4 square degrees of sky in each exposure. The telescope 
is capable of viewing to a limiting magnitude of R ~ 20. The 
QUEST team devoted twelve 60 s exposures over several nights 
for each trigger in both the winter and autumn periods, with a 
target rate of 1 trigger per week. 

ROTSE III (lAkerlofetal.l2003h is a collection of four robotic 
telescopes spread around the world, each with a 0.45 m aperture 
and 3.4 square degree FOV. No filters are used, so the spectral 
response is that of the CCDs, spanning roughly 400 to 900 nm. 
The equivalent R band limiting magnitude is about 17 in a 20 s 
exposure. The ROTSE team arranged for a series of thirty images 
for the first night, and several images on following nights, for 
each autumn run trigger, with a target rate of 1 trigger per week. 

SkyMapper ( iKeller et al.ll2007h is a survey telescope located 
at Siding Spring observatory in Australia. The mosaic camera 
covers 5.7 square degrees of sky in each field, and is mounted 
on a 1.35 m telescope with a collecting area equivalent to an 
unobscured 1.01 m aperture. It is designed to reach a limiting 
magnitude ^ ~ 21 (>7 sigma) in a 110 s exposure. SkyMapper 
accepted triggers in the autumn run with a target rate of 1 per 
week, with several fields collected for each trigger. 

TAROT jKlotz et al.ll2009al) operates two robotic 25 cm tele- 
scopes, one at La Silla in Chile and one in Calern, France. Like 
the ROTSE III system, each TAROT instrument has a 3.4 square 
degree FOV. A 180 second image with TAROT in ideal condi- 
tions has a limiting R magnitude of 17.5. During the winter run, 
TAROT observed a single field during one night for each trig- 
ger. In the autumn run, the field selected for each trigger was 
observed over several nights. TAROT accepted triggers with a 
target rate of 1 per week. 



Liverpool j 



V- 




;SKy I 



I SkyMapper 



Fig. 1. A map showing the approximate positions of telescopes 
that participated in the project. The Swift satellite observatory is 
noted at an arbitrary location. The image is adapted from a blank 
world map placed in the public domain by P. Dlouhy. 



Zadko Telescope (ICoward et al.ll2010l) is a 1 m telescope lo- 
cated in Western Australia. The current CCD imager observes 
fields of 0.15 square degrees down to magnitude ~ 20 in the 
R band for a typical 180 s exposure. For each accepted trigger 
in the autumn run, Zadko repeatedly observed the five galaxies 
considered most likely to host the source over several nights. The 
target trigger rate for Zadko was one trigger per we ek. 

The Liverpool telescope (Steele et af] l2004l) is a 2 m 
robotic telescope situated at the Observatorio del Roque de Los 
Muchachos on La Palma. For this project the RATCam instru- 
ment, with a 21 square arcminute field of view, was used. This 
instrumentation allows a five minute exposure to reach magni- 
tude r' - 21. This project was awarded 8 hours of target-of- 
opportunity time, which was split into 8 observations of 1 hour 
each, with a target rate of 1 trigger per week. 

3.2.2. Radio and X-ray Instruments 

LOFAR dFender et al.l 120061: Ide Vos et alj 120091: IStappers et all 
[2011') is a dipole array radio telescope based in the Netherlands 
but with stations across Europe. The array is sensitive to fre- 
quencies in the range of 30 to 80 MHz and 1 10 to 240 MHz, and 
can observe multiple simultaneous beams, each with a FWHM 
varying with frequency up to a maximum of around 23°. During 
the autumn run, LOFAR accepted triggers at a target rate of 1 
per week and followed up each with a four-hour observation in 
its higher frequency band, providing a ~25 square degree field 
of view. 

Although not used in the prompt s earch during the sc ience 
run, the Expanded Very Large Array (iPerlev et alj|201 ll) was 
used to follow up a few triggers after the run with latencies of 
3 and 5 weeks. 

The Swift satellite dGehrels et al.ll2004l) carries three instru- 
ments, each in different bands. Swift granted several target of 
opportunity observations with two of these, the X-ray Telescope 
(XRT) and UV/Optical Telescope (UVOT), for the winter and 
autumn observing periods. The XRT is an imaging instrument 
with a 0.15 square degree FOV, sensitive to fluxes around 10"'^ 
ergs/cm^/s in the 0.5-10 keV band. A few fields were imaged for 
each trigger that Swift accepted. 



4. Trigger Selection 

The online analysis process which produced GW candidate trig- 
gers to be sent to telescopes is outlined in Fig. |2] After data and 
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Fig. 2. A simplified flowchart of flie online analysis with approx- 
imate time requirements for each stage. Data and information 
on data quality were generated at the Hanford, Livingston, and 
Virgo interferometers (HI, LI, and VI) and copied to central- 
ized computer centers. The online event trigger generators pro- 
duced coincident triggers which were written into the GraCEDb 
archive. The LUMIN and GEM algorithms selected statistically 
significant triggers from the archive and chose pointing loca- 
tions. Significant triggers generated alerts, and were validated 
manually. If no obvious problem was found, the trigger's esti- 
mated coordinates were sent to telescopes for potential follow- 
up. 



information on data quality were copied from the interferome- 
ter sites to computing centers, three different data analysis algo- 
rithms identified triggers and determined probability skymaps. 
The process of downselecting this large collection of triggers to 
the few event candidates that received EM follow-up is described 
in this section. 

After event candidates were placed in a central archive, addi- 
tional software used the locations of nearby galaxies and Milky 
Way globular clusters to select likely source positions (Sect.|5]l. 
Triggers were manually vetted, then the selected targets were 
passed to partner observatories which imaged the sky in an at- 
tempt to find an associated EM transient. Studies demonstrating 
the performance of this pipeline on simulated GWs are presented 
in Sect.m 

4.1. Trigger Generation 

Sending GW triggers to observatories with less than an hour la- 
tency represents a major shift from past LIGO/Virgo data anal- 
yses, which were reported outside the collaboration at soonest 
several months after the data collection. Reconstructing source 
positions requires combining the data streams from the LIGO- 
Virgo network using either fully coherent analysis or a coinci- 
dence analysis of single-detector trigger times. A key step in la- 
tency reduction was the rapid data replication process, in which 
data from all three GW observatory sites were copied to several 
computing centers within a minute of collection. 

For the EM follow-up program, three independent GW de- 
tection algorithms (trigger generators), ran promptly as data 
became available, generating candidate triggers with latencies 
between five and eight minutes. Omega Pipeline and coherent 
WaveBu rst (cWB), which are both described in lAbadie et al.l 
(1201 Obi) , searched for transients (bursts) with only loose as- 
sumptions regarding waveform morphology. The Multi-Band 
Template Analysis (MBTA) (Marion 2004), searched for sig- 
nals from coalescing compact binaries. Triggers were ranked by 
their "detection statistic", a figure of merit for each analysis, 
known as Q, t], and Pcombined, respectively. The statistics 77 for 



cWB and Pcombined for MBTA are related to the amplitude SNR 
of the signal across all interferometers while Q is related to the 
Bayesian likelihood of a GW signal being present. Triggers with 
a detection statistic above a nominal threshold, and occurring in 
times where all three detectors were operating normally, were 
recorded in the Gravitational-wave Candidate Event Database 
(GraCEDb). 

The trigger generators also produced likelihood maps over 
the sky (skymaps), indicating the location from which the signal 
was most likely to have originated. A brief introduction to each 
trigger generator is presented in Sects. 143711 - 14.1.31 

4.1.1. Coherent WaveBurst 

Coherent Wave Burst has been used i n pre vious searche s for GW 
bursts, such as lAbbott etaP (l2009bl) and lAbadie et all (1201 Obi) . 
The algorithm performs a time-frequency analysis of data in the 
wavelet domain. It coherently combines data from all detectors 
to reconstruct the two GW polarization waveforms li+{t) and 
hx{t) and the source coordinates on the sky. A statistic is con- 
structed from the coherent terms of the maximum UkeUhood ra- 
tio functional dPlanagan & Hugheslll998l:lKhmenko et al.ll2005h 
for each possible sky location, and is used to rank each lo- 
cation in a grid that covers the sky (skymap). A detailed de- 
scription of the likelihood analysis, the sky localization statistic 
and th e performance of the cWB algorithm is published else- 
where ( lKlimenkoetalJl2011h . 

The search was run in two configurations which differ in 
their assumptions about the GW signal. The "unconstrained" 
search places minimal assumptions on the GW waveform, while 
the "linear" search ass umes the signal is dom inated by a single 
GW polarization state ( Kli menko et al.ll2bl ih . While the uncon- 
strained search is more general, and is the configuration that was 
used in previous burst analyses, the linear search has been shown 
to better estimate source positions for some classes of signals. 
For the online analysis, the two searches were run in parallel. 

4.1.2. Omega Pipeline 

In the Omega Pipeline search ( lAbadie et alj 1201 Obh . triggers 
are first identified by performing a matched filter search with 
a bank of basis waveforms which are approximately (co)sine- 
Gaussians. The search assumes that a GW signal can be decom- 
posed into a small number of these basis waveforms, and so is 
most sensitive to signals with a small time-frequency volume. 
Coincidence criteria are then applied, requiring a trigger with 
consistent frequency in another interferometer within a physi- 
cally consistent time window. A coherent Bayesian position re- 
construction code ( Searle et al. 2008, 2009J is then applied to 
remaining candidates. The code performs Bayesian marginaliza- 
tion over all parameters (time of arrival, amplitude and polariza- 
tion) other than direction. This results in a skymap providing the 
probability that a signal arrived at any time, with any amplitude 
and polarization, as a function of direction. Further marginaliza- 
tion is performed over this entire probability skymap to arrive at 
a single number, the estimated probability that a signal arrived 
from any direction. The Q. statistic is constructed from this num- 
ber and other trigger properties. 

4.1.3. MBTA 

The Multi-Band Template Analysis (MBTA) is a low-latency 
implementation of the matched filter search that is typically used 
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to sea rch for compact binary inspirals (lMarionll2004t iBuskulid 
l2oToh . In contrast to burst searches which do not assume any 
particular waveform morphology, MBTA specifically targets the 
waveforms expected from NS-NS, NS-BH and BH-BH inspi- 
rals. In this way it provides complementary coverage to the burst 
searches described above. 

The search uses templates computed from a second order 
post-Newtonian approximation for the phase evolution of the 
signal, with component masses in the range 1-34 Mq and a total 
mass of < 35 Mq. However, triggers generated from templates 
with both component masses larger than the plausible limit of the 
NS mass — conservatively taken to be 3.5 Mq for this check — 
were not considered for EM follow-up, since the optical emis- 
sion is thought to be associated with the merger of two neutron 
stars or with the disruption of a neutron star by a stellar-mass 
black hole. 

Triggers from each interferometer are clustered and used to 
search for coincidence among the individual detectors. To gen- 
erate a candidate event for follow-up, triggers with consistent 
physical parameters must be present in all three LIGO/Virgo in- 
terferometers. For each triple coincident trigger, the sky location 
was estimated using the time delay between detector sites and 
the am plitude of the signal measured in each detector ( iFairhurstI 
l2009l) . Before the observing period, a set of simulated gravita- 
tional wave signals was used to measure the distribution of er- 
rors in recovering the time delays and signal amplitudes. The 
sky localization algorithm then uses these distributions to assign 
probabilities to each pixel on the sky. 



4.2. Estimating False Alarm Rates 

The primary quantity used to decide whether an event should be 
considered a candidate for follow-ups was its FAR, the average 
rate at which noise fluctuations create events with an equal or 
greater value of the detection statistic. For the winter run, a FAR 
of less than 1 event per day of livetime was required to send an 
imaging request to the ground-based telescopes, with a higher 
threshold for Swift. For the autumn run, the FAR threshold was 
0.25 events per day of livetime for most telescopes, with stricter 
requirements for sending triggers to Palomar Transient Factory 
and Swift. Livetime is here defined as time all three interferom- 
eters were simultaneously collecting usable science data. 

As in p revious all-sky burst searches, e.g. lAbbott et al.l 
(I2009bl) and lAbadie et al.1 (1201 Obi) , the FAR for the two burst 
pipelines was evaluated using the time-shift method. In this 
method, artificial time shifts, between one second and a few hun- 
dred seconds, are applied to the strain series of one or more in- 
terferometers, and the shifted data streams are analyzed with the 
regular coherent search algorithm. The shifted data provide an 
estimate of the background noise trigger rate without any true 
coincident gravitational wave signals. During the online anal- 
ysis, at least 100 time shifts were continuously evaluated with 
latencies between 10 minutes and several hours. The FAR of 
each event candidate was evaluated with the most recent avail- 
able time shifts. 

The MBTA pipeline evaluated the FAR analytically based on 
single interferometer trigger rates, rather than using time shifts. 
This is computationally simpler than the burst method. It is valid 
since MBTA is a coincident rather than a coherent analysis, and 
allows the FAR to be evaluated with data from the minutes im- 
mediately preceding the trigger time ( lMarionll2004l) . 



4.3. Online Data Quality 

A number of common occurrences may make a stretch of inter- 
ferometer data unsuitable for sensitive GW searches. Examples 
include times of large seismic disturbance, non-standard inter- 
ferometer configurations, and temporary saturations of various 
photodiodes in the interferometer sensing and control system. 
To mark such times, monitor programs analyze auxiliary data to 
produce lists of abnormal time segments with low latency. When 
a trigger was identified, it was automatically checked against 
these lists; triggers which occurred in stretches of unacceptable 
data were automatically rejected. During this search, all three 
GW detectors were simultaneously collecting science quality 
data for roughly 45% of the time. 

4.4. Manual Event Validation 

In addition to automated checks on data quality, significant trig- 
gers were manually vetted. Trigger alerts were broadcast to col- 
laboration members via e-mail, text message, a website, and in 
the interferometer control rooms as audio alarms. For each alert, 
a low-latency pipeline expert conferred with personnel at each of 
the three observatory sites to validate the event. Pipeline experts 
and scientists monitoring data on-site provided 24/7 coverage in 
8 hour shifts. Assigned personnel confirmed the automated data 
quality results, checked plots for obvious abnormalities, and ver- 
ified that there were no known disturbances at any of the three 
observatory sites. 

The intention of manual event validation was to veto spuri- 
ous events caused by known non-GW mechanisms that have not 
been caught by low-latency data quality cuts, not to remove ev- 
ery non-GW trigger. In fact, at current sensitivities, most or all 
of the triggers are unlikely to represent true astrophysical events. 
The trade-off for this additional check on the quality of the trig- 
gers was added latency (usually 10 to 20 minutes) between trig- 
ger identification and reporting to the EM observatories. It is 
possible that as the search matures in the Advanced LIGO/Virgo 
era the validation process can be fully automated. 



5. Choosing Fields to Observe 

The uncertainty associated with GW posi tion estirnates, ex- 
pected to be several tens of square degrees ( lFairhurstll2009l) . is 
large compared to the FOV of most astronomical instruments. 
Moreover, the likely sky regions calculated from interferometer 
data may be irregularly shaped, or even contain several disjoint 
regions. It is impractical to image these entire regions given a 
limited amount of observing time for a given instrument. There 
is thus a need to carefully prioritize fields, or tiles, of an instru- 
ment to optimize the likelihood of imaging the true gravitational 
wave source. 

The LUMIN software package was created to gather GW 
triggers from the three trigger generators, and use the skymaps 
and locations of known galaxies to select fields for each opti- 
cal or radio instrument to observe. In addition, LUMIN includes 
tools that were used to facilitate trigger validation (Sect. 14.41 ) and 
communication with robotic telescopes. Fields for observation 
with the Swift XRT and UVOT were selected with slightly dif- 
ferent criteria by a separate software package, the Gravitational 
to Electro-Magnetic Processor (GEM). During the testing pro- 
cess, GEM also applied the tiling criteria for optical telescopes 
to simulated skymaps, and so provided an important consistency 
check between LUMIN and GEM. 
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5.1. Galaxy Catalog 

The G ravitational Wave Galaxy Catalog (GWGC) dWhite et alJ 
was created to help this and future searches quickly iden- 
tify nearby galaxies. 

The catalog contains up-to-date information compiled from 
the literature on sky position, distance, blue magnitude, ma- 
jor and minor diameters, position angle and galaxy type for 
53,225 galaxies rangi ng out to 100 Mpc, as well as 150 Milky 
Way globular clusters. [White et all (1201 ll) compared the catalog 
with an expected blue light distribution derived from SDSS data 
and concluded that the GWGC is nearly complete out to ~40 
Mpc. The catalog improves on the issue of multiple entries for 
the same galaxy suffered by previous catalogs by creating the 
GWGC from a subset of 4 large catalogs, each of which lists 
a uniqu e Principal Galaxy Catalogue (PGC) number for every 
galaxy (iPaturel et al.lll989l) . The catalogs used were: an upd ated 
version of the Tully Nearby Galaxies Catalog (Tullv'l987h, the 
Catalog of N eighboring Galaxi es (Karachentsev et al. 2 0041). the 
V8k catalog dTullv et al.ll2009l) . and HyperLEDA dPaturel et alJ 
I2003h . Also included is a Ust of 150 known Milky Way globular 
clusters These are all available freely online, but a 

local, homogeneous list is essential for rapid follow-up purposes. 

5.2. Weighting and Tiling Algorithm 

To make use of the galaxy catalog, and choose tiles for each GW 
trigger, similar algorithms have been implemented in the GEM 
and LUMIN software packages. 

The position information from the trigger generators (see 
Sect. 14.1b is encoded in skymaps that assign a likelihood to each 
0.4° X 0.4° pixel in a grid covering the sky. In practice, only 
the 1000 most likely pixels are retained, limiting the sky area to 
roughly 160 square degrees. The search volume is further lim- 
ited by keeping only objects in the catalog with an estimated dis- 
tance of less than 50 Mpc, as the current sensitivity of the GW 
detectors makes it unlikely that binaries containing a neutron star 
would be detectable beyond this distance. Approximately 8% of 
the pixels in an average skymap contain a local galaxy or globu- 
lar cluster listed in the GWGC catalog. 

For burst triggers, the tiling algorithms treat the luminosity 
of each galaxy or globular cluster as a prior for its likelihood 
to host a GW emitting event. The blue light luminosity is used 
as a proxy for star formation, indicating the presence of mas- 
sive stars that may be GW burst progenitors themselves and may 
evolve into compact binaries that eventually merge. In addition, 
weak sources of GWs are assumed to be more numerous than 
strong sources, so that a closer galaxy should contain more de- 
tectable sources than a more distant galaxy of the same mass 
dNut tall & Sutton 2010J. This leads to assigning the following 
likelihood to each pixel: 



P oc 



Y Mil 



(1) 



where L is the likelihood based only on the GW data, and M 
and D are the blue light luminosity (a rough proxy for mass) and 
distance of the associated galaxy or globular cluster The sum is 
over all the objects associated with a particular pixel (which will 
be or 1 galaxy for the majority of pixels). Extended nearby 
sources which have a major axis larger than the pixel size have 
their mass divided evenly over each pixel falling within the el- 
lipse of the disk defined by their major and minor axes. Once 
this calculation is performed for each pixel, the entire skymap is 
renormalized to a total likelihood equal to unity. 



Unlike the burst algorithms, MBTA assumes the GW source 
is a merging binary, and estimates some of the source's physical 
parameters for each trigger This allows the galaxy catalog to be 
applied in a slightly different way. Each interferometer measures 
a quantity known as effective distance 



£>eft = D 



1 + cos- L 



+ fI cos^ t 



-1/2 



(2) 



where D is the actual distance to the source, l is the inclination 
angle between the direction to the observer and the angular mo- 
mentum vector of the binary, and F+ and Fx are the antenna re- 
sponse functions of the particular interferometer The important 
feature of the effective distance is that it is always greater than or 
equal to the true distance to the source. For each MBTA trigger 
the galaxy catalog is then only considered out to the smallest ef- 
fective distance measured for that trigger, with a maximum pos- 
sible effective distance of 50 Mpc. After the catalog is downse- 
lected in this way, each pixel is weighted by the fraction of the 
catalog's total mass contained in that pixel, i.e. 



p^y Mf'^'L, 



(3) 



with the sum over all galaxies associated with the pixel, and 
2i Mf"'= = 1 for a sum over the downselected catalog. 

These procedures require a pixel's coordinates to be consis- 
tent with a known galaxy's location to be targeted by telescopes. 
However, in the case that the skymap does not intersect with 
any galaxies in the catalog, the likelihood from the GW skymap 
alone is used as each pixel's likelihood {P = L). In practice, this 
is a very rare occurrence and only happens in the case of a very 
well-localized skymap. 

The actual pointing coordinates requested for each telescope 
are selected to maximize the total contained P summed over pix- 
els within the FOV. If multiple pointings are allowed with the 
same instrument, additional tiles with the next highest ranking 
are then selected. The tile selection process is illustrated in Fig. 

s 



5.3. Galaxy Targeting for Small-Field Instruments 

The logic used for selecting pointings for the Swift satellite was 
similar to that of ground-based telescopes, except that, because 
the narrower Swift FOV required greater precision, care was 
taken to ensure the target galaxies were within the selected field. 
The coordinates supplied to Swift for follow-up were those of 
the matched galaxy itself in cases where there was only a single 
galaxy in a pixel, but the center of the 0.4°x0.4° pixel in cases 
where the central coordinates of an extended source were outside 
the pixel or there were multiple galaxies in the pixel. Since fewer 
follow-ups were allowed using Swift than with other scopes, a 
minimum requirement was placed on the statistic P contained 
within the pixels selected for X-ray observation. 

Zadko and Liverpool Telescope also have relatively narrow 
fields. For these telescopes, no attempt was made to capture mul- 
tiple galaxies in a single field. Instead, the weighting scheme in 
Eqn.[T]was applied to each galaxy rather than each pixel, and the 
center coordinates of the top ranked galaxies were passed to the 
observatories. 
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Fig. 3. The weighting and tiling process for a simulated signal reconstructed by cWB. The skymap is shown in the left panel with 
the highest likelihood regions in red, and lower ranked pixels in blue, along with galaxy locations marked as black circles. The right 
panel shows the location and approximate size of the three chosen QUEST tiles, along with the locations of pixels that are retained 
after weighting by the galaxy catalog. The injection location is caught by the southernmost tile, and is marked with an asterisk. 



6. Observing Strategy 

6.1. Communication 

After an event candidate passed manual inspection, a script was 
launched to pass the GPS time and selected field center loca- 
tions to the QUEST, ROTSE 111, SkyMapper, TAROT, Zadko, 
Liverpool Telescope, and LOFAR observatories. During the au- 
tumn run, a total of five such alerts were sent. During the (ear- 
lier) winter run, 8 event candidates were passed to the TAROT 
and QUEST observatories. The number of field locations passed 
to each telescope for each GW event candidate are listed as 
the "Tiles per Trigger" in Table |2] During the autumn run, in 
cases where the fields selected for a particular instrument were 
unobservable due to daylight or latitude, no alert was sent to 
the observatory. In most cases, alerts were sent via a direct 
socket connection from a LIGO computing center at Caltech 
with IP mask protection. Alerts to ROTSE 111, SkyMapper, 
TAROT, and Zadko used the format of GCN notices. Alerts to 
LOFAR and the Liverpool Telescope used the VOEvent format 
(IWilliams & Seamanl2006h . For QUEST, the GPS time and field 
positions were posted as ASCII tables to a password protected 
web site which was regularly polled by the QUEST scheduler. 

The Palomar Transient Factory received field locations and 
GPS times using the VOEvent format via a socket connection, 
but with a more restrictive FAR threshold than the other optical 
telescopes, and so triggers were only sent to PTE if the on-call 
team executed a separate script. Alerts to Swift also required ex- 
tra action by the on-call team, who entered field coordinates in 
an online form. The Pi of the Sky prototype telescope was en- 
gaged through automated e-mails and manual checks of a pass- 
word protected web page. 

6.2. Telescope Response 

The wide variety of telescopes involved in the search led to a 
diversity of observing strategies, with each partnering group ap- 
plying a different cadence. By design, most of the telescopes in 
the network were robotic, and could respond to alerts without 



human intervention. In a few cases this allowed response times 
of less than a minute after an alert was sent, though response 
times of a few hours were more typical due to wait time for tar- 
gets to be overhead. 



During the winter run, QUEST responded to three triggers, 
making 2 exposures of each field on the night of the request. 
TAROT responded to one winter run trigger, taking six images 
on the night of the request. Swift also responded to one trigger 
in the winter run, taking one exposure of each field following the 
request, and then a second set of exposures on a later date to be 
used as reference images. 



For most observatories in the summer run, the observing plan 
called for capturing a first image of the selected fields as rapidly 
as possible, with follow-up observations every night or every 
other night out to five days after the trigger time. For the op- 
tical observatories, any night's observation included 2 or more 
exposures for each field, to help eliminate asteroids, CCD ar- 
tifacts, and other contaminants from the data set. In addition, 
some fields were imaged at later times, up to a month after the 
trigger time, to provide reference images, or possibly to capture 
a light curve with a late brightening time. TAROT, Zadko, PTE, 
QUEST, and Pi of the Sky all followed this recipe. ROTSE exe- 
cuted a more aggressive observing plan, collecting a set of 30 im- 
ages in rapid succession on the first night, and then sets of eight 
images on each of 15 nights following the trigger with intervals 
of two days on average. As in the winter run. Swift made one 
exposure of each field following the trigger, and then collected 
a reference image after a lag of several weeks. The Liverpool 
Telescope devoted roughly one hour of observation upon receiv- 
ing a trigger, and then collected reference images a few weeks 
after the trigger time. The LOFAR response was not automated. 
A telescope operator made a single, four hour observation one 
to four days after delivery of a trigger. SkyMapper also required 
manual intervention to respond to a trigger, and so responded on 
a best effort basis. 
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Fig. 4. Plots of typical uncertainty region sizes for the Omega (top) and unconstrained cWB (bottom) pipelines, as a function of 
GW strain amplitude at Earth, for various waveform types. The "searched area" is the area of the skymap with a likelihood value 
greater than the likelihood value at the true source location before the galaxy catalog is used to further limit the search region. The 
solid line with symbols represents the median (50%) performance, while the upper and lower dashed lines show the 75% and 25% 
quartile values. Near the detection threshold {h^^s ~ 10"^' Hz"'^^) , uncertainty regions are typically between 10 and 100 square 
degrees. The Omega pipeline performs poorly on white noise bursts but exceptionally well on sine-Gaussians because it is designed 
to identify signals that are well-localized in frequency space. 



7. Performance Study 

7.1. Simulated Waveform Injections 

An ensemble of simulated GW signals was generated to mea- 
sure the effectiveness of the reconstruction and follow-up pro- 
cedures. For the Omega and cWB burst pipelines, these "soft- 
ware injections" were a mix of ad hoc sine-Gaussian, Gaussian, 
and white noise burst waveforms similar in type and distri- 
bution to those used in previous LIGO/V irgo all-sky analyses 
dAbbott et al.ll2009bt IXbadie et al.ll2010"bl) . While these wave- 
forms are not based on specific astrophysical models, they do 
a good job of characterizing detector response for signals in 
specific frequency ranges (sine-Gaussians) and broadband sig- 
nals (white-noise bursts). ForMBTA (see Sect. 14. 1.3b . injections 
were drawn from NS-NS and NS-BH inspiral waveforms with 
a range of parameters. To emulate a realistic spatial distribu- 
tion, each injection was calculated with a source distance and 
direction inside a randomly selected galaxy from the GWGC 
and the simulated GW amplitudes were weighted to be inversely 
proportional to distance. Only galaxies within 50 Mpc were in- 
cluded in the simulation, with weighting factors applied so that 
the probability of originating from each galaxy was proportional 



to its blue light luminosity. The simulation set and the analysis 
used the same catalog, so the results presented in Figures |6]- [8] 
make the assumption that the blue light luminosity distribution 
of galaxies in the GWGC is a good tracer of GW sources in the 
local universe. Signals were superimposed on real LIGO- Virgo 
gravitational wave data taken between August and December 
2009. 

While performance studies in this paper were done using 
software injections, a relatively small number of tests in which 
a signal was physically put into the interferometer via actuators 
("hardware injections") were also performed, providing an addi- 
tional cross-check. 

7.2. Testing Results 

Because the skymap likelihood regions are often irregularly 
shaped, the size of the uncertainty region is characterized by the 
"searched area", defined as the angular area of the skymap with 
likelihood greater than the likelihood at the true source location. 
The median searched area as a function of signal strength is plot- 
ted for both cWB and Omega Pipeline in Fig. |4] Here, gravita- 
tional wave amplitudes are expressed in terms of their root-sum- 
squared amplitude: 
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Fig. 5. Plots of uncertainty region sizes for the MBTA pipeline 
as a function of combined SNR (pcombined)- The solid line with 
symbols represents the median (50%) performance, while the 
upper and lower dashed lines show the 75% and 25% quartile 
values. The expected detection threshold is around Pcombined ~ 
12. 
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(4) 



where h+ (t) and hx (0 are the plus- and cross -polarization strain 
functions of the wave. Since h is a dimensionless quantity, /2,-ss 
is given in units of Hz For this data, signals near the de- 
tection threshold would have h^ss ~ 10'^ ' Hz~'^^, roughly cor - 
responding to the cWB statistic 77-5 dAbadie et al ll2010bh . 
These signals were typically localized with median search areas 
of several tens of square degrees. The coherent position recon- 
struction algorithms are "tuned" to localize these near-threshold 
signals as accurately as possible; as a result, some of the plots 
reveal a degradation in algorithm performance for very loud sig- 
nals. Median searched area is shown for MBTA in Fig. |5] as a 
function of the combined SNR of the signal: 



Pcombined — 'sjPlil "^^Ll Py 



(5) 



where Php /'lp P\\ '■^^ single detector SNRs seen in the 
Hanford, Livingston and Virgo instruments, respectively. 

The simulated GW signals described above were also used 
to test the tiling software in order to determine the success rate 
for imaging the correct sky location with realistic detector noise, 
reconstructed skymap shapes, and telescope FOVs. The LUMIN 
software package was used to determine pointings for ground- 
based telescopes and GEM was used for Swift. 

Some of the results of this simulation can be seen in Fig. 
|6] The results are plotted as a function of the ranking statistic 
used by each pipeline. On the y-axis, the "Fraction of triggers 
imaged" represents the fraction of triggers with the given detec- 
tion statistic that have the correct image location included within 
the selected tiles. Given a GW trigger, the success rate plotted 
in Fig. |6] estimates the odds of choosing the right sky position. 
In this figure, note that the "whole skymap" is limited to 160 
square degrees, and so does not always include the true source 
location. The thresholds for initiating follow-ups varied with the 
condition of the interferometers, but was typically around 3.0 
for Q, 77 = 3.5 for cWB, and Pcombined = 10 for MBTA. The 



complex behavior of the Omega efficiency curve is related to the 
use of a hybrid detection statistic which utilizes different meth- 
ods depending on SNR range. Clearly, events with SNR near the 
threshold for triggering follow-ups, the most likely scenario for 
the first detections, are the most difficult to localize. 

Example efficiency curves for the burst simulation are shown 
in Fig.|7] The efficiency for each marker on the plot is calculated 
as the fraction of signals for which the injected location was 
successfully imaged, for an /irss range centered on the marker. 
Specifically, we require that: 

1. The trigger's ranking statistic is higher than the threshold, 
which is chosen to enforce a FAR of about 1 GW trigger per 
day of livetime. 

2. The true source location is included in one of the chosen 
tiles. Five tiles are allowed for Swift, three tiles for the 
QUEST camera, and one tile for all other telescopes. 

Note that efficiencies in this figure do not reach unity even for 
loud events primarily due to the difficulty of correctly localizing 
GWs in some regions of the sky where the antenna response of 
one or more interferometers is poor. 

The efficiencies produced with these criteria are upper Hmits 
on what would be detected in a real search. They assume that 
the EM transient is very bright, and will always be detected if 
the cotTect sky location is imaged. The quoted efficiencies do not 
account for the fact that some chosen tiles will not be observed 
due to restrictions from weather, instrument availability, proxim- 
ity to the Sun or Moon, or the application of a manual veto. The 
exact behavior of the efficiency curves will vary depending on 
the morphologies of the simulation waveforms selected. Finally, 
the chosen GW trigger FAR of one event per day presumes the 
false alarm fraction from the EM transient classification pipeline 
will be low enough to make a coincidence significant. Many of 
these additional complications and their associated impact on ef- 
ficiency are described in Metzger & Berger (20111). 

Nevertheless, these curves provide a measure of the potential 
for joint EM/GW searches. If the number of incidental EM tran- 
sients in the observed fields can be understood and controlled, 
then the addition of EM data can effectively increase the search 
sensitivity to very weak GW signals. For occasional strong GW 
signals, the plots suggest that only a few pointings of a telescope 
are enough to image the true location with better than 50% effi- 
ciency. 

7.3. Calibration Uncertainty 

Uncertainty in the calibration of GW detectors ( lAbadie et al.l 
2010a; Accadia et al. 2011) may impact the ability to correctly 
choose the right fields to observe with EM instruments. To esti- 
mate the potential detriment to pointing, we generated a second 
set of simulated burst signals, with each signal including some 
level of miscalibration corresponding to realistic calibration er- 
rors. Before being added to detector noise, each astrophysical 
signal was scaled in amplitude by a factor between 0.85 and 
1.15, and shifted in time by between -150 and 150 yus. The exact 
amplitude and time "jitter" were randomly selected from flat dis- 
tributions for each signal entering each detector. The bounds of 
the distribution of values for the timing and amplitude jitter were 
chosen to match preliminary estimates for the LIGO and Virgo 
calibration error budgets around 150 Hz for the 2009-2010 run. 
Well above this frequency, the actual timing errors are likely less 
than this model; the simulation is conservative in this sense. 

Some of the results of this simulation, with the cWB algo- 
rithm, may be seen in Fig. [8] The success rate is shown for the 
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Fig. 6. Success rates for the tile selection process based on unconstrained cWB (left). Omega (right), and MBTA (bottom) skymaps. 
An injection recovered with the detection statistic shown on the horizontal axis is considered a success if the correct source location 
is included in one of the chosen tiles. Typical thresholds for follow-up are Q=3.0, ?7=3.5, and Pcombined=10. Each tile is 1.85° x 1.85°, 
the FOV of both the ROTSE and TAROT telescopes. Statistical uncertainties are small with respect to the markers. 



entire pipeline, assuming one pointing of a 1.85° x 1.85° FOV, 
three pointings of the QUEST FOV, and five pointings of a Swift 
FOV. The curves are shown both with and without the effects 
of calibration uncertainty. For the low SNR signals that are the 
most likely for first detections, 77 < 10, the efficiency is within 
a few percent with and without the calibration uncertainty. This 
is expected, since at low signal to noise ratio, timing uncertainty 
from de tector noise is la rger than timing uncertainty due to cali- 
bration (lFairhurstll2009h . However, for louder signals, the ability 
to correctly choose the right sky location is seen to be modestly 
impacted by the accuracy of the calibration. 

8. Summary 

Mergers of compact binary systems containing neutron stars, as 
well as some other energetic astrophysical events, are expected 
to emit observable transients in both the gravitational wave and 
electromagnetic channels. Observing populations of joint signals 
would likely reveal many details of the GW sources, and could 
even constrain cosmological models. 

During 2009 and 2010, the LIGO and Virgo collaborations 
partnered with a large, heterogeneous group of EM observato- 
ries to jointly seek transient signals. X-ray, optical, and radio 
observatories collected follow-up observations to GW triggers 



that were delivered with ~30 minutes of latency. Analysis of the 
multi-instrument data set is currently in progress, and the results 
of the search for jointly observed transients will be published at 
a later date. 

A Monte Carlo study of the GW data analysis algorithms 
used in the low latency pipeline demonstrated the ability of the 
LIGO/Virgo network to localize transient GW events on the sky. 
Localization ability depends strongly on the SNR of the GW 
signal; lower SNR signals are more difficult to localize, but are 
also the more likely scenario for the first detections. Signals with 
SNR near the detection threshold were localized with median 
sky areas between 10 and 100 square degrees. After limiting 
the search to known galaxies and Milky Way globular clusters 
within the detection range of the GW observatories, the simu- 
lation shows that the correct location of signals detected near 
threshold can be imaged with 30-50% success with three fields 
of size 1.85° x 1.85°, for instance. Moreover, the ability to im- 
age the source position is seen to be only marginally impacted 
by realistic levels of calibration uncertainty. 

This search establishes a baseline for low-latency analysis 
with the next-generation GW detectors Advanced LIGO and 
Advanced Virgo. Installation of these second-generation detec- 
tors is already in progress, with observations expected to begin 
around 2015. Developing a low-latency response to GW triggers 
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Fig. 7. Fractional success as a function of strain at Earth for combinations of the Omega and cWB burst pipelines. Success rates 
assume 5 pointings for each event for Swift (left) and 3 for QUEST (right). Fractional success is end-to-end from triggering by 
pipeline to successful pointing, with a threshold for follow-up approximating a FAR of one per day. The "Logical Or" curve counts 
a success if either linear cWB or Omega correctly localized the event, effectively doubling the allowed number of tiles. Some curves 
show degraded performance for very loud signals because the algorithms are tuned to optimize performance for weaker events close 
to the detection threshold. Statistical uncertainties are small with respect to the markers. 




Fig. 8. Fractional success rate for simulated gravitational wave signals with (dashed) and without (solid) including calibration 
uncertainties. QUEST results assume 3 pointings. Swift results assume 5 pointings and TAROT/ROTSE results assume 1 pointing. 



represents the first steps toward solving the many logistical and 
technical challenges that must be overcome to collect prompt, 
multiwavelength, EM observations of GW source progenitors. 
The integration of GW and EM observatories is likely to con- 
tinue to develop over the next few years as the scientific com- 
munity prepares to utilize the many opportunities promised by 
the impending global network of advanced GW detectors. 
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Name 


nana 


FOV (square degrees) 


Aperture (m) 


Exposure Time (s) 


Limiting Ma 


Palomar Transient Factory 


Optical 


/.J 


l.Z 


oU 


zU.J 


or toe i)ky 


Uptical 


4UU 


f\ mo 
U.U/Z 


lU 


1 1 


QUEST 


Optical 


9.4 


1 


60 


20 


T •} i' x 1 " C I ; TTT 

K(J1 sSb 111 


Optical 


1 A 

3.4 


0.45 


zU 


1 T C 


Sky Mapper 


Optical 


J. 1 


1 .jj 




Z i 


TAROT 


Optical 


3.4 


0.25 


180 


17.5 


Zadko Telescope 


Optical 


0.15 


1 


180 


20 


Liverpool Telescope 


Optical 


0.0058 


2 


3600 


21 


LOFAR 


Radio 


-25 


N/A 


14400 


N/A 


Swift 


X-ray 


0.15 


N/A 


200-5000 


N/A 


Swift 


UV, Optical 


0.078 


0.3 


200-5000 


24 



ignitude 



Notes. Some characteristics of instruments involved in the search. The shown limiting magnitudes are estimates, assuming favorable observing 
conditions. 

Table 2. Partner Instrument Follow-up Information 



Name 


Run 


Tiles per Trigger 


Target Alerts Per Week 


Triggers Imaged 


Palomar Transient Factory 


Autumn 


10 


1/3 


1 


Pi of the Sky 


Autumn 


1 


1 


1 


QUEST 


Both 


3 


1 


5 


ROTSE 111 


Autumn 


1 


1 


5 


SkyMapper 


Autumn 


~9 


1 


3 


TAROT 


Both 


1 


1 


3 


Zadko Telescope 


Autumn 


5 


1 


2 


Liverpool Telescope 


Autumn 


1 


1 


1 


LOFAR 


Autumn 


1 


1 


2 


Swift 


Both 


5 


1/4 


2 



Notes. Follow-up information for instruments involved in the search. Each instrument participated in either the autumn run, or both the winter 
and autumn runs. The column marked "Tiles per Trigger" shows how many different field locations the instrument attempted to observe for each 
accepted GW event candidate. The 'Target Alerts Per Week" column shows that alerts were sent to PTF and Swift at a lower rate than the other 
observatories. The final column shows the number of GW event candidates for which each instrument collected data. 
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